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SU^IMARY 

The coupling of the combined source-vortex distribution of Green’s poten- 
tial flow function with contemporary numerical techniques is shown to provide 
accurate, efficient, and stable solutions to subsonic inviscid analysis and 
design problemc for multi-element airfoils* The analysis problem is solved by 
direct calculation of the surface singularity distribution required to satisfy 
the flow tangency boundary condition. The design or inverse problem is solved 
by an iteration process. In this process, the geometry and the associated 
pressure distribution are iterated until the pressure distribution most nearly 
corresponding to the prescribed design distribution is obtained. Typically, 
five iteration cycles are required for convergence. A description of the 
analysis and design method is presented, along with supporting examples. 

INTRODUCTION 

The surface panel method philosophy for solving arbitrary incompressible 
potential flow problems involves the mating of classical potential theory with 
contemporary numerical techniques. Classical theory is used to reduce an 
arbitrary flow problem to a surface integral equation relating boundary condi- 
tions to an unknown singularity distribution (Reference l),The contemporary 
numerical technique", are then used to calculate an approximate solution to the 
integral equation (References 2-14). 

/J.1 properly formulated surface panel methods are exact in the sense that 
the difference between the approximate numerical solution and the exact solu- 
tion to the integral equation can be made arbitrarily small at the expense of 
increasing the number of computations. This does not imply that all panel 
methods are equally successful. Indeed, vast differences exist with respect 
to prediction accuracy versus computational effort, reliability, simplicity, 
and applicability to an inverse solution mode for design problems. 

The major distinguishing characteristics of panel methods are depicted in 
Figure 1. For the special case of two-dimensional flow, nearly all the pos- 
sible combinations of Figure 1 have been formulated and tested at McDonnell 
Aircraft Company (MCAIR). The most successful of the tested formuLations was 
selected as the foundation for the MCAIR Multi-Element Airfoil Analysis and 
Design Computer Program, herein after designated the MCAIR Airfoil Program. 

The analysis (direct) mode calculates the velocity distribution of an arbi- 
trary airfi il geometry, whereas the design (inverse) mode iterates to generate 
the geometry most nearly corresponding to a prescribed surface velocity distri- 
bution. Unusually rapid and consistent convergence is obtained because the 
inverse algorithm includes all the first order terms in the relationship be- 
tween arbitrary geometry and velocity perturbations. 
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This paper presents ^background theory, the formulation, and representa-^ 
tive numerical solutions for the MCAIR Airfoil Program. 


NOMENCLATURE 


Airfoil lift coefficient 
c Pressure coefficient [1 - (^ )^] 

n Unit normal vector to a boundary 

V Flow velocity 

Y Vortex density 

y Doublet density 

a Source Density 

SUBSCRIPTS 

E External (fluid) side of a boundary 

I Internal side of a boundary 

N Normal component 

T Tangential component 

Free stream conditions 

SURFACE SINGULARITY THEORY 

Any three-dimensional, incompressible, potential flow field can be con- 
sidered to be induced by a suitable distribution of source and doublet singu- 
larity densities on flow boundary surfaces (Reference !)• It is usually con- 
venient to treat the flow field as the sum of a uniform free stream Voo plus a 
disturbance potential field associated with the presence of the body (Figure 
2). Then the total velocity V at any field point can be expressed as 

V = + ^<1. (1) 


where (fi is the potential of the disturbance field. 

The value of (Ji at an arbitrary field point P can be expressed in the 
following form: 
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( 7 )} ds 


where ds is a differential area element at arbitrary surface point Q of the 
body; r is the distance to point P; n is the distance measured along an axis 
nor^l to the body surface at Q, positive outward; and 0 and p are functions 
of location on the body surface, i.e., functions of Q. a and v are respective- 
ly the source and doublet densities* 

For two-dimensional flow, it is generally more convenient to use vortex 
singularities in place of doublets* As is proved in Reference 4, a surface 
doublet distribution of density p can be replaced by an equivalent vortex 
distribution where the vortex density vector y satisfies the following equa- 
tion at each surface point: 


Y - n X Pj 

n is the local unit normal vector pointing ii:to the flow field. 

It is noteworthy that there is no limit '.j the number of different solu- 
t on source - vortex distributions correspondxng to any given flow field. The 
theoretical distinction between the different distributions is best lllus- 
trated by examining the imaginary flow field internal to the boundaries (Flg- 
Consider tvo— dimensional flow. The discontinuity across any surface 
sheet of sources and Vortices -an be expressed as follows: 


VNe - VNj - o 

% - = Y (5) 

where o and y are the local source and vortex densities. For solid body bound- 
ary conditions (Vne = 0), Equation (4) indicates that a vortex-only solution 
will correspond to Vui = 0 at every internal boundary point. The unique in- 
ternal flow field generated by zero normal velocity boundary conditions is, 
of course, stagnation. Then Equation (5) indicates that VTp = v everywhere 
on the extent^ suriace. 

A particularly useful combii^ed source-vortex distribution corresponds to 
the uniform internal flow field Vi - V„. For this case, Equations (4) and (5) 
imply that the source and vortex densities are equal to the external pertur- 
batlon velocity components, i*e* , 

\ - Vn„ + o 

\ - Vt„ + Y 

Subscript E has been omitted for brevity. The above combined source-vortex 
distribution is equivalent to the application of Green’s third Identity (Refer- 
ence 15) to the perturbation potential (}i. A schematic of Equations ( 6 ) and 
(7) is presented in Figure 4. v / uu 


( 6 ) 

(7) 
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An illustration of the nature of three different singularity distributions 
corresponding to the same potential field is presented in Figure 5. The theo- 
retical singularity distributions for a source solution, vortex solution, and 
the combined Green s identity source-vortex solution are shorn for the flow 
aroun an infinite circular cylinder with flow tangency boundary conditions. 

It IS interesting that the combined source-vortex distribution is a fifty per- 
cent blend of the source only and vortex only solutions, Green's identity 
t^ically provides a source distribution more mild than the source only solu- 
tion and a vortex distribution more mild than the voxtex only solution. 

^e combined source-vortex distribution of Green's identity is especially 
suitable for application to a numerical panel method. The equality between 
singularity densities and disturbance flow velocity means that regardless of 
either boundary conditions or geometry, the singularity magnitudes cannot be- 
come disproportionately large. This contrasts sharply with source only solu- 
tions, tor which the source density can increase without bounds as body thick- 
ness approaches zero even though the velocity remains finite everywhere. The 
Significance is that the mild source-vortex distributions associated 
® identity eliminate excessive velocity gradients between boundary 
condition control points, thereby reducing the possibility of leakage. The 
relationships between disturbance flow velocity and singularity density holds 
regardless of boundary conditions, be they Neumann, Dirlchlet, or mixed. 

These relationships can always be used to eliminate half the unknown singular- 
ity densities a priori, leaving no more effective unknowns that a source only 
or vortex only approach. ^ 


TWO-DIMENSIONAL SOLUTION FORMULATION 

formulation is based on the combined source-vortex distriba- 
tion of Green s identity (Equations 6 and 7). Low order panel modeling is 
emp oye to the effect that source gradient and surface curvature corrections 
on each panel are ignored. However, the use of internal potential boundary 

application of the velocity - singularity strength equality 
results in close to higher order prediction accuracy for most practical^geo- 
mctrxc shapes. The advantage of the low order modeling is the simplicity 
inherent in establishing inverse capability. 

The geometry of each airfoil element is simulated by a closed polygon, 
^ere the polygon segment end points are assumed to lie on the actual airfoil. 
The midpoint of each segment (panel) is selected as the boundary condition 


For Neumann prescribed normal velocity boundary conditions, 
density distribution is established a priori from Equation (6), 
that the source density is uniform on each panel. 


the source 
It is assumed 


The gist of the solution to analysis problems is to determine the appro- 
priate vortex distribution. The vortex density is assumed to vary linearly 
on each panel and to be continuous at panel end points, with one exception. 

If the geometry has any slope discontinuities such as a sharp trailing edge, 
the vortex density is allowed to be discontinuous at the corresponding panel 
end point. This is consistent with exact theoretical solutions, for which 


224 




f- 


®® discontinuities in the vortex distribution. The 
mgnitude of the discontinuity is an additional unknow which is determined 

introduction of an internal boundary condition control point near 

local, panel lSgth“‘'“‘'® nom:aally selected as 1% of the 

<*irectly imposing prescribed normal velocity boundary condi- 
tioM at the control points, a theoretically equivalent approach firL applied 

cLf S:id t consistent with^L 

city field yi = V„, uniform internal perturbation potential is prescribed at 

lnt^!i ^ accomplished by specifying that the analytical line 

^ ^ perturbation velocity component along an internal path connecting 

adjacent control points be zero. connecting 

The circulation of each airfoil element is controlled by either of two 

Tr'^ctlZ d‘S.«lT Tf ‘“e^gtS^rbe 

Jh! 7..^^ Alternately, a Kutta condition can be applied in which 

the velocity normal to the trailing edge bisector is set eq^«l to zero It a 
location approximately 2% local segment length downstream. The nominal 27 
value has been select.J Ucause it cons^..tently provides Iifrcoe?SLt pre- 

Of catherall et al. (Reference 16) for typical panel models. 

Imposition of the potential boundary conditions and one circulation con- 
P« airfoil element establishes a system of lineareq^tJon^ 
the same number of unknown vortex densities as equations. This number is equal 
to the sum of the total number of panels and total number of shaTcoroers! 
lutlon to the system renders the complete set of singularities known. 

(7) a^d surface velocity is calculated from Equation 

LLulIf, calculated by 



( 8 ) 


Force and moment integration is performed under the 
trol point pressure applies uniformly to each panel. 


assumption that the con- 


izflMnn solutions are generated in accordance with the iterative-linear- 
zation philosophy developed in Reference 17. For the elimentstobe designed 
a multi-element airfoil systenv-the following steps are Involved: 


( 1 ) 


The user prescribes a design pressure or 
the surfaces of the various elements. 


velocity distribution 


around 
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(2) The user prescribes a starting geometry to initialize the calcula- 
tions and the location of one point per element to be fixed in space, 
such as the trailing edge. 

(3) The program solves the direct problem for the geometry, 

in order to determine the change in velocity distribution required 
to achieve the prescribed values. 

(4) The program calculates the rate of change of surface velocity with 
respect to an arbitrary change in surface angle distribution. Each 
element perimeter remains fixed. If the tangential conponent of 
velocity at the control point of the ith panel is designated Vti and 
if the surface angle of the jth panel is designated 0j, then the 
array is calculated where 

3Vt^ 

^Ij " 'W' 

(5) The change in surface angle distribution is calculated in accordance 
with the prescribed velocities and the following first order expres- 
sion; 

= E A0j) 

(6) The geometry is corrected by the program and steps (3)- (5) are 
repeated as a series of iteraclon cycles. 


The most difficult and Important step in formulating the Inverse capability is 
to generate the matrix A^j. It is noted that all terms were incorporated in 
deriving the partial derivative, including singularity strength changes and 
the displacement of panels j+1, j+2, etc. corresponding to the surface angle 
change d©j. The corresponding singularity strength changes are obtained by a 
first order expansion to the boundary condition equation. In order to obtain 
numerical stability in the Inverse process, the velocities are prescribed not 
only at the panel midpoints but at endpoints as well. Then the solution in- 
volves minimizing the mean square error between the prescribed and calculated 
distributions. 


EXAMPLE SOLUTIONS 

The following examples demonstrate the numerical behavior of the present 
MCAIR Airfoil Program as compared to alternate approaches. 


The benefits of using a combined source-vortex distribution versus source- 
only is improved accuracy for (1) thin highly loaded airfoils and for 
(2) sharp concave corners. Examples of the frrner are presented in Reference 
14 and are not repeated here. The latter is typical of wing-fuselage inter- 
sections and is represented by the wedge-cylinder example of Figure 6. Both 
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Douglas Neumann source method (Reference 2) were 
^ geometry using identical panel modeling (244 panels) The 

re^il sharp^cLLHorne;, which 

in thr^p,?^TJT internal potential boundary conditions is reflected 

Figure 7 Karman-Trefftz airfoil of 

^ ^ ^ calculated panel method distributions the 

vartlnJ^f^M Catharall-Sells solution was prescribed. Con- 

gradient and panel curvature effects Fnr r. 4 ^ 4 -u^ . source 

cUrwarconilterbr^Se ^’'°t singularity densUy^S^u^facf velo- 

Dofenf-iai K j ^ j^® most accurate approach, namely. Equation (7) for 
florM^gew” conditions and the sunmation of influence coefficients for 

The third example demonstrates prediction accuracv of ^ha mpatd -^-u j 

r^Jese'StJT^ 

corresponding to the trailing edge h/sector’ 

r.„n I**® Objective of the first Inverse example is to design a circular 
cylinder by the MCAIR method using a nearly flat nlata 

fFieure tho ^ i3t plate for the starting geometry 

and the co;verg^ solution geoLtrrof"FiguS°9"wL‘^ob^aiSefa^ 

radius of endpoints are within a maximum distance of 0 002 

iui resLct derivatives of vei;cJty 

ohtfli 1 ^ surface angle change were necessary but not sufficient for 

comM"^® convergence about the periphery of this exan.ple. usi of Sid 

combined source-vortex singularites is also a factor. To illustraL the 

SZe%rflTa:?d aingularltiua uarr^ed^u 

±3S5r c^rjSg :^ge ^isii£;nc!££HF 
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analysis solution using 70 panels total requires approximately 4 seconds com- 
puting time. For the corresponding inverse solution, approximately 20 seconds 
is required per iteration cycle, 

CONCLUDING REMARKS 

The use of the source-vortex distribution of Green's identity coupled 
with Internal perturbation potential boundary conditions provides a simple, 
accurate, reliable, and efficient procedure for solving airfoil incompressible 
potential flow problems. For a wide range of geometric shapes, low order, 
paneling modeling generates prediction accuracy usually associated with higher 
order solutions, The procedure is especially suitable for application to in- 
verse design problems. 
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TABLE 1.- MODIFICATIONS OF THE 


Subroutine 

No 

change 

MAIN 


POINT 

X 

SLOPE 

X 

TRANS 


DISTP 

X 

FTLUD 


DIR 

X 

LSQ 


PROOT 

X 

MAIN1 


READIT 


GEOM 


ROTRAN 


ASLOT 


NORMAL 


MAIN 2 


CHEN 


MATRIX 


POTLF 


CAMBER 


SMOOTH 

X 

VOVBT 


THICK 


COMPR 


STAG 


MAINS 


LOAD 


LAMNA 


BLTRAN 


TURBL 


TURB 


OERIV 


START 


CONFBL 


CONF5 


CONF7 


CONFB 


DLIM 

X . 
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3-. 
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SINGULARITY TYPE 


SOLID BODY BOUNDARY CONDITIONS 


• SOURCE • ZERO NORMAL VELOCITY 

• VORTEX (DOUBLET) • CONSTANT STREAM FUNCTION 



Figure 1.- Diagram of surface singularity panel methods for 
Incompressible* potential flow. 



Figure 2 


Body Immersed In an unbounded flow field 
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Figure 4.- Green's Identity velocity-singularity relationship. 
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Figure 9.- Circular-cylinder Inverse solution. MCAIR method. 


Sterling Gtometry 




Figure 10.- Circular-cylinder inverse solution. 
Vortex-only method. 



